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TESTS OF EVALUATED BERYLLIUM (n,2n) CROSS SECTIONS BY ANALYSIS OF 
1.4-eV FLUX IN WATER FOR AN AMERICIUM-BERYLLIUM SOURCE 
ENCLOSED IN BERYLLIUM SPHERES 

by Daniel Fieno, Roger L. Alexander, and C. Hubbard Ford 
Lewis Research Center 

SUMMARY 

Measurements in water of the neutron flux distributions at indium resonance energy 

(1.4 eV), with beryllium (Be) shells enclosing an americium-beryllium (Am-Be) neutron 

source, were made to quantitatively test two evaluated sets of (n,2n) cross sections for 
9 

Be. Multigroup Sn calculations of the experimental flux distributions have been made 

Q 

using cross sections from these two evaluations of the decay modes of the Be(n,2n) reac- 
tion. The first evaluation is based on the data files of the GAM-II program which use the 
same cross section evaluation employed for the current ENDF/B library for Be. The 
second evaluation is based on Perkins' detailed studies of the reactions which contribute 

Q 

to the Be(n,2n) cross section. 

This work required an accurate evaluation of the spectrum of the neutron source. An 
uncollided source spectrum was obtained which, when used in an S n calculation, gave 
agreement with the measured flux distribution for the source in water. The neutron leak- 
age spectrum computed for the source in air using this uncollided spectrum was also re- 
quired to match a measured leakage spectrum. 

Q 

Tests of the evaluated Be(n,2n) cross sections are made by analysis of measure- 
ments in water of the flux distributions for the Am-Be source enclosed in Be spheres. 

The calculations using cross sections for Perkins' Be(n,2n) evaluation gave good agree- 
ment with the experiments. Calculations with Be(n,2n) cross sections obtained from the 
data files of the GAM-II code gave poor agreement with these experiments. 

For the Am-Be source in the 1. 9- centimeter-thick Be shell, the computed enhance- 
ments at 1.4-eV are 9.4 and 9. 5 percent for calculations using Perkins' Be(n,2n) data 
and GAM-II Be data, respectively. The experimental enhancement is measured to be 
7.3±2.0 percent. For the Am-Be source in the 7. 0- centimeter-thick Be shell, the com- 
puted enhancements are 26.3 and 27.6 percent for calculations using Perkins' Be(n,2n) 
data and GAM-II Be data, respectively. The experimental enhancement is measured to 
be 22. 8±2. 5 percent. 

The calculation of the shape of the spatial distribution of 1.4-eV energy neutrons and 
values of enhancement using Perkins' Be(n,2n) data are significantly better than using the 
GAM-II Be(n,2n) data. Since current ENDF/B cross sections are based on the same 
evaluation as used to formulate the GAM-II cross sections, it is suggested that Perkins’ 
evaluation be used to update cross sections for Be. 


INTRODUCTION 


Analysis of reactors using beryllium (Be) or beryllium oxide as either a moderator 
or a reflector requires accurate calculation of the so-called fast multiplication effect due 
to (n, 2n) reactions. Precise criticality calculations have been hampered by an incom- 
plete knowledge of the decay modes of the Be(n, 2n) reaction and large discrepancies in 
both the fast effect and age still exist (ref. 1). Since beryllium is of considerable im- 
portance to nuclear reactors, it is essential to assess by means of additional experi- 
ments the relative applicability of various models for the Be(n, 2n) reaction. 

In the present experiments measurements of the neutron flux distribution at indium 
resonance energy were made in water for an americium -beryllium (Am-Be) neutron 
source enclosed by beryllium shells. Multigroup transport theory calculations of these 
experimental flux distributions are compared using two evaluations of the Be(n, 2n) 
cross sections. Thus, the calculated as compared to the measured flux distributions 
provide quantitative indications of the relative effects of various cross section evalua- 
tions of the decay modes of the Be(n, 2n) reaction. 

Two Be(n, 2n) cross section evaluations are considered for this study. The first 
evaluation is based on a widely used cross section pregram, GAM-II (ref. 2), which 
grossly approximates the Be(n, 2n) reaction mechanism in computing the multigroup 
neutron transfer cross sections. The GAM-II cross sections are based on an evalua- 
tion of Joanou and Stevens (ref. 3) which also serves as a basis for the present ENDF/B 
cross section library for Be (ref. 4). The second evaluation is based on the work by 
Perkins described in references 5 to 7 which presents detailed evaluations of the 
Be(n, 2n) cross sections and provides calculations of the angular and energy distributions 

Q 

for neutrons from the various decay modes of the Be(n, 2n) reaction. This work by 
Perkins provides a more comprehensive evaluation of the multigroup transfer cross 
sections for beryllium. 

In the first part of the present study an uncollided or virgin spectrum was con- 
structed for the Am-Be source by means of S n transport calculations for two experi- 
ments. In the first experiment a leakage or collided spectrum for the source in air was 
measured using an NE-213 proton recoil scintillator. In the second experiment the in- 
dium resonance energy flux distribution was measured for the source in an effectively 
infinite water medium. By an iterative procedure, a virgin spectrum was obtained 
which correctly calculated this flux distribution in water, as well as reproduced the 
measured source leakage spectrum in air. This virgin spectrum was required not to 
exceed the estimated limits of error of published small source spectra measured by 
several independent techniques. 

The second part of this work employs the Am-Be source, enclosed by spherical 
shells of Be, in an infinite medium of water to measure the neutron flux distributions 
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at indium resonance energy. Multigroup S n transport calculations are compared to 
corresponding experimental flux distributions and to source enhancements due to Be. 


DESCRIPTION OF EXPERIMENTS AND REDUCTION OF DATA 

Flux distributions at indium resonance energy were measured in water for three 
different spherical configurations: for the Am-Be source alone and for the source at 
the center of each of two beryllium spheres. 


Source 

12 

The spherical Am-Be neutron source contains 54 curies (2x10 disintegrations/sec) 
24 1 

of Am (16. 7 g) in a net volume of 84. 3 cubic centimeters. There are 4. 0±0. 4 grams 

241 

of beryllium for each gram of Am, thus giving a total of 66. 8±6. 7 grams of Be in 

O 

the source. The output of the source was (1. 30±0. 08)xl0 neutrons per second in Jan- 
uary 1967. 

The outer surface of the source is spherical (6. 096-cm outside diameter) and in- 
cludes two stainless steel containment shells, one 0. 140 centimeter thick and the other 
0. 152 centimeter thick. The dimension of .the spherical cavity for the source material 
is 5. 462 centimeters inside diameter. The Am-Be material was pressed in the inner 
two hemispheres and a spring was placed between the two before welding. 


Tank 

For all experiments except the bare source in water, the indium resonance energy 
flux distributions were measured in a cylindrical tank of demineralized water that was 
2. 44 meters outside diameter by 2. 44 meters high. The tank used for the bare source 
in water experiment was 1. 83 meters outside diameter by 1. 52 meters high. Both 
tanks provided an effectively infinite water medium for these configurations. 


Beryllium Spheres 

The two spheres used were the following: 

(1) 10. 16-centimeter outside diameter Be sphere (1. 9-cm shell thickness) 

(2) 20. 32-centimeter outside diameter Be sphere (7. 0-cm shell thickness) 

All spheres were made of hemispheres. The beryllium hemispheres were made with a 


3 



6. 35-centimeter inside diameter central space to accommodate the source. To center 
the source and also displace any water, an aluminum spherical shell 0. 127 centimeter 
thick was placed around the source. 


Foils and Holder 

To measure the neutron flux distributions at indium (In) resonance energy, indium 
foils 4. 99 centimeters outside diameter by 0. 0127 centimeter thick were used. The 
large diameter foils were required so that measurements could be made as far as 
80 centimeters from the bare source in water. These foils were covered all around by 

cadmium 0. 089 centimeter thick corresponding to a cutoff energy of 0. 6 eV. The ac- 

1 1 

tivity of the irradiated foils is then proportional to the 54-minute half-life in ac- 

115 

tivity. The 4. 4-hour half-life In activity caused by inelastic scattering of neutrons 
above about 1 MeV was small because the high energy flux decreases rapidly as a func- 
tion of distance in water and because irradiations close to the source were of short dur- 
ation of the order of 30 minutes. Measurements as far as approximately 100 centime- 
ters from the source were made with these foils in aluminum covers. 

A lucite holder was used to position the foils. Figure 1 shows the foil holder po- 
sitioned over a sphere. The three legs were notched to locate the cadmium covered 
foils. The notches were spaced every 1. 27 centimeters for the first 11. 43 centimeters 
and every 2. 54 centimeters from there on out to nearly 150 centimeters. No two foils 
were ever placed closer than 7. 62 centimeters apart during any one irradiation. 


Counters 

Each foil was counted on both sides in several 2n beta counters. Counter plateaus 

were obtained periodically using irradiated indium foils and carbon 14 (^C). The use 

of several counters provided a check of counter stabilities during the counting periods. 

Counter efficiencies were not measured. To reduce background counts, one of the 

hemispherical counters of radius 2. 86 centimeters was equipped with an anticoincidence 

scintillation region which surrounded all but the bottom of the counter. A minimum- 

ionizing muon will dissipate at least 3 MeV in this scintillator. An approximately 1-MeV 

threshold was obtained by setting the scintillator voltage at the point where cobalt 60 
fi 0 

( Co) gammas began to be detected. Thus, most background counts caused by muons 
are subtracted. This particular counter background was reduced to about 1.7 counts per 
minute from about 23 counts per minute. Foil backgrounds were counted in each counter 
prior to each irradiation. This background measurement becomes important for those 
foils located in the water far from the source. 
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Data Reduction 


Counting data were reduced using the following equation: 


C 


s, m 



T 


C + L - B 




( 1 ) 


where 

B foil background 

C raw counts 

O 

L dead time counting loss, (t^/t c )C 

t count time, min 

t^ dead time, min 

tj irradiation time, min 

t w time from end of irradiation to start of count, min 
t half-life/ln 2 (mean life, min) 

C g is the measured saturated foil count rate as opposed to a saturated activity 
since counter efficiencies were not measured. The average C_ computed for the 

S y Ixl 

foil front sides was averaged with the average C e „ for the foil back sides. This 
average C g m is taken to be proportional to the indium resonance energy flux (ref. 8). 

Foil size corrections were made to account for the geometric variation of fluxes 
over the various segments of each foil. Different portions of these large foils were 
at different distances from the center of symmetry. The data were corrected using the 
following equation (ref. 9): 


C fi (r n ) = C t 


's' r 0' = ^s, m^O* - ^T 


l d C s,mM 
dr 


where 


C (r n ) measured saturated foil count rate at r n 

dC s m (r 0 )/dr approximation to flux derivative at Tq 

r f 
r 0 


( 2 ) 


foil radius, cm 

distance from source center to center of foil, cm 
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IIIIIIHIIIIIII 


TABLE I. - EXPERIMENTAL INDIUM RESONANCE ENERGY ACTIVITY FOR Am-Be 
NEUTRON SOURCE IN INFINITE WATER MEDIA 


[Data normalized using equation (4b).] 



Source alone 


10. 16-cm o. d. Be sphere around 

20. 32-cm o. d. Be sphere around 





source a 



source 3 " 


Radius, 

Corrected 

Error 

Radius, 

Corrected 

Error 

Radius, 

Corrected 

Error 

cm 

counts, 


cm 

counts, 


cm 

counts, 



c s 



C s 



C s 


3. 46 

0. 1470E-03 

0. 59E-06 

5. 50 

0. 1287E-03 

0. 23E-06 

10. 57 

0. 9801E-04 

0. 52E-06 

4.73 

. 1377E-03 

. 11E-06 

6.67 

. 1150E-03 

. 50E-06 

11. 84 

. 7813E-04 

. 38E-06 

6.00 

. 1152E-03 

. 10E-06 

8.04 

. 9449E-04 

. 16E-05 

13. 11 

. 5651E-04 

. 28E-06 

7.27 

. 9015E-04 

. 92E-07 

9.31 

.7215E-04 

. 79E-06 

14.38 

. 3877E-04 

. 24E-06 

8. 54 

. 6730E-04 

. 52E-07 

10. 57 

. 5154E-04 

. 22E-06 

15.65 

. 2585E-04 

. 37E-06 

9. 81 

. 4982E-04 

. 62E-07 

11. 87 

. 3674E-04 

. 45E-06 

16. 92 

. 1730E-04 

. 17E-06 

11.08 

. 356 IE-04 

. 39E-07 

13. 12 

. 2672E-04 

. 16E-06 

18. 19 

. 1169E-04 

. 38E-07 

12. 35 

. 2670E-04 

. 36E-07 

14.39 

. 1907E-04 

. 67E-07 

19.46 

. 8212E-05 

. 46E-07 

14. 89 

. 1463E-04 

. 45E-07 

15.66 

. 1411E-04 

. 12E-06 

20.73 

. 5797E-05 

. 45E-07 

17. 43 

. 8336E-05 

. 18E-07 

16.94 

. 1033E-04 

. 57E-07 

22.00 

. 4182E-05 

. 26E-07 

19. 97 

. 5046E-05 

. 17E-07 

19. 48 

. 5747E-05 

. 10E-06 

24. 54 

. 2349E-05 

. 20E-07 

22. 51 

. 3 10 IE-05 

. 10E-07 

22.00 

. 3529E-05 

. 41E-07 

27.08 

. 1383E-05 

. 16E-07 

25.05 

. 1895E-05 

. 28E-07 

27.09 

. 1299E-05 

. 35E-07 

29.62 

. 845 IE-06 

. 56E-08 

30. 13 

. 7755E-06 

. 38E-08 

32. 17 

. 5259E-06 

. 76E-08 

32. 16 

. 533 IE-06 

. 15E-07 

32.67 

. 5052E-06 

. 26E-08 

37. 25 

. 2352E-06 

. 49E-08 

34. 70 

. 3401E-06 

. 25E-08 

35. 21 

. 3446E-06 

. 27E-08 

39.79 

. 1540E-06 

. 31E-08 

37. 24 

. 2236E-06 

. 33E-08 

37. 75 

. 23 14E-06 

. 34E-08 

44. 87 

. 7517E-07 

. 60E-08 

39. 78 

. 1476E-06 

. 38E-08 

40. 29 

. 1507E-06 

. 55E-08 

49. 95 

. 3541E-07 

. 15E-08 

42. 32 

. 9876E-07 

. 94E-09 

42. 83 

. 10 13E-06 

. 17E-08 

52. 49 

| . 2247E-07 

. 17E-08 

44. 86 

. 6606E-07 

. 10E-08 

45. 37 

. 6685E-07 

. 11E-08 

57. 57 

! . 1165E-07 

. 96E-09 

47. 40 

. 4601E-07 

. 11E-08 

47. 82 

. 5086E-07 

. 11E-08 

65. 19 

. 3732E-08 

. 28E-09 

49. 94 

. 3068E-07 

. 16E-08 

50.45 

. 3322E-07 

. 12E-08 

67. 73 

. 2672E-08 

. 73E-10 

52. 48 

, . 2151E-07 

. 72E-09 

52. 99 

. 2301E-07 

. 64E-09 

70. 27 

. 1826E-08 

. 37E-09 

55.02 

. 1477E-07 

. 29E-09 

55. 53 

. 16 19E-07 

. 48E-09 

72. 81 

. 9449E-09 

. 26E-09 

57. 56 

. 1039E-07 

. 37E-09 

58.07 

. 1096E-07 

. 42E-09 




60. 10 

. 7258E-08 

. 39E-09 

60.61 

. 7731E-08 

. 27E-09 




62. 64 

. 4477E-08 

. 76E-09 

63. 15 

. 6108E-08 

. 83E-09 




65. 18 

. 3475E-08 

. 45E-09 

65. 69 

. 4349E-08 

. 37E-09 




67. 72 

. 2538E-08 

. 44E-09 

68.23 

. 3273E-08 

. 21E-09 




70. 26 

. 1458E-08 

. 24E-09 

70. 76 

. 2203E-08 

. 32E-09 







73. 27 

. 1291E-08 

. 35E-09 







75. 84 

. 9255E-09 

. 23E-09 







80. 92 

. 5192E-09 

. 24E-09 









a Source is centered in Be spheres by 0. 127-cm-thick aluminum spherical shell. 
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Counting Statistics 


Counting statistics were calculated in several ways. Theoretical errors, that is, 
the square root of the number of counts, were obtained for each count. Also for each 
counter a standard deviation of the saturated count rate was calculated. Composite er- 
rors were calculated over all counters for both the theoretical and standard errors. 
These composite errors were then compared to a standard deviation of all the C g *s 
for a particular foil as though the results from all counters were part of the same dis- 
tribution. Counter plateau shifts or noise problems could be seen easily in this way. 
Composite errors were obtained in all cases with the C g from front and back sides 
of the foil combined to give an average value C g m> 

Finally, when measurements at a given radius were repeated, the standard devia- 
tion of the C_ __ and a composite error were again computed. The larger error is 

S y IXX 

reported here. The distribution data and the errors are listed in table I for the three 
experiments with the Am-Be neutron source. 


CALCULATIONAL METHODS AND CROSS SECTIONS 

Multigroup calculations using the S n formulation of transport theory were used to 
obtain the indium resonance energy flux distributions for all of the experiments consid- 
ered. The formulation of this method is given in an appendix of reference 10. The cal- 
culations were all done in spherical geometry with an effectively infinite water region, 
the outer boundary radius being greater than 110 centimeters. The Am-Be neutron 
source and its stainless steel container were explicitly included in the calculations both 
as to dimensions and composition. The calculations were for a fixed source using a 
virgin spectrum and a uniform distribution throughout the source region. 

Thirty energy groups were used to describe the spectrum for the various configura- 
tions considered. The details of this group structure are presented in table II. For 
convenience, group 28 which extends in energy from 1. 125 to 1. 445 eV, was chosen as 
the group representative of the indium resonance activation energy. 

Two computer programs were used to provide the group averaged microscopic 
cross sections for all of the isotopes used in the calculations. These programs are 
GAM -II (ref. 2) and GATHER-II (ref. 11). The GAM-II program also provided one of 
the evaluations for the P Q down scattering transfer cross sections for the Be(n, 2n) 
reaction. 

GAM-II is a Bg multigroup code for the calculation of fast neutron spectra. For 
this study, the code was used to integrate the Bg equations over 99 energy groups 
covering the energy range from 14. 92 MeV to 0. 414 eV. The cross sections for the 
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TABLE II. - 30-GROUP ENERGY AND 


illllllllllllllllllllllll 


LETHARGY BOUNDARIES 


Group 

Upper energy 
boundary, 
eV a 

Lethargy 

boundary 

Lethargy 

width 

1 

L 105X10 7 

-0. 1 

0. 1 

2 

1. 000X10 7 

0 



3 

9. 048X10 6 

. 1 



4 

8. 187X10 6 

.2 



5 

7. 408x10® 

.3 



6 

6. 703X10® 

.4 



7 

6. 065X10® 

. 5 



8 

5. 488X10® 

.6 



9 

4. 966X10® 

.7 



10 

4. 493X10® 

.8 



11 

4. 066X10® 

.9 


2 

12 

3. 329X10® 

1. 1 


2 

13 

2. 725X10® 

1.3 


3 

14 

2. 019X10 6 

1.6 


4 

15 

1. 353X10 6 

2.0 


5 

16 

8. 209X10 5 

1 2. 5 


7 

17 

4. 076X10 5 

3.2 


7 

18 

2. 024X10 5 

3. 9 


85 

19 

8. 652X10 4 

4.75 

1 

50 

20 

1. 931X10 4 

6.25 

1 

75 

21 

3. 355X10 3 

8.0 

1.25 

22 

961. 1 

9.25 

1. 

0 

23 

353.6 

10.25 



24 

130. 1 

11.25 



25 

47. 85 

12.25 



26 

17.60 

13. 25 

1. 

25 

27 

5.043 

14. 50 

1. 

25 

28 

1. 445 

15.75 

. 

25 

29 

1. 125 

16.0 

1. 

0 

30 

.414 

17.0 

-- 

-- 


a Thermal group (group 30) lower energy 
boundary is 0. 001 eV. 


29 fast groups are obtained by averaging over the 99-group spectrum obtained for a 
given material composition. Down scattering transfer cross sections for elastic scat- 
tering through the P« order are included for all isotopes. 

GATHER-II was used to obtain thermal, or 30 in group, cross sections. This pro- 
gram is a Pj code for the calculation of thermal spectra over 101 energy points cover- 
ing the energy range from 2. 38 to 0. 001 eV. All of the cross sections obtained using 
the GATHER-II program are based on a temperature of 293 K. 
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The Pq down scattering transfer matrix for the Be(n, 2n) reaction in the GAM-II 
data files is based on the evaluation of Joanou and Stevens (ref. 3). This evaluation 
gives a total cross section for the Be(n, 2n) reaction as well as partial cross sections for 
the excitation of the 2. 43- and 1. 75-MeV levels. These cross sections are shown in 
figure 2(a) for comparison with the evaluated Be(n, 2n) cross section data of Perkins 
(ref. 5) which is shown in figure 2(b). Both evaluations of the total Be(n, 2n) cross sec- 
tions are in agreement. However, the GAM-II evaluation of the decay modes of the 
Be(n, 2n) reaction are in disagreement with the evaluation due to Perkins. 

The largest contribution to the Be(n, 2n) reaction in Perkins’ evaluation comes from 

Q 

the excitation of the 2. 43-MeV level in Be*. The majority (91. 7 percent) of these re- 

O 

actions are assigned as leaving Be* in a broad 2. 90-MeV excited state with an effective 

O 

excitation energy in Be of 0. 6 MeV. The remaining reactions go to the ground state 

8 9 

of Be. Smaller cross sections for the excitation of levels in Be* at 6. 76 and 9. 1 MeV 

are also assigned by Perkins. 

The remaining cross section which is a large fraction of the Be(n, 2n) cross section 

at lower energies proceeds by multibody breakup. Two multibody models for energy 

and momentum transfer are ordinarily considered. One is the cluster model. In this 

9 8 

model, after inelastic scattering of a neutron with Be, the residual Be* nucleus and 
the second neutron move together, finally separating. This model is indistinguishable 

Q 

from the excited state reaction in which the Be* nucleus is left in the 1. 70-MeV ex- 

O 

cited state. In this case, the excitation energy of the Be* nucleus shared with the sec- 
ond neutron is small (35 keV). The second model partitions the energy in proportion to 
the momentum phase volume. The reactions of interest are the three-body breakup 
Be(n, 2n) Be* and the four-body breakup Be(n,2n)2o!. This report assumes that the 

O 

multibody breakup to the Be ground state proceeds by the three-body phase space model 
rather than the cluster model. 

The evaluation of the reported Be(n, 2n) data by Perkins (ref. 5) was used to com- 
pute the Be(n, 2n) Pq down scattering transfer matrix using the AGN- SIGMA code 
(ref. 7). In the present report, this evaluation of the Be(n, 2n) data is referred to sim- 
ply as Perkins’ Be data. The Pq down scattering cross sections based on GAM-II are 
radically different from Perkins' Be data. Figure 3 shows a comparison of the Pq 
down scattering data for several neutron energies using the GAM-II evaluation' and Per- 
kins' evaluation. For an incident neutron energy of 9. 5 MeV the GAM-II Pq down 
scattering data as shown in figure 3(a) extends down to 3. 7 MeV while Perkins' Be data 
extends below 0. 1 MeV. Similar results obtained for incident neutron energies of 5. 8 
and 3. 5 MeV are shown in figures 3(b) and (c). 

For the transport calculations the number of energy groups, the elastic scattering 
order, and the type and order of the quadrature set used in the S n calculations are im- 
portant. Mynatt, Greene, and Engle (ref. 12) found that elastic scattering treated to 
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TABLE HI. - GAUSS- LEGENDRE 


QUADRATURE SET FOR S 16 
TRANSPORT CALCULATIONS 


IN SPHERICAL GEOMETRY 


Direction 

Direction cosine 

Weight 

number 



1 

-1.0 

0 

2 

-.9894009 

.0135762 

3 

-.9445750 

.0311268 

4 

-.8656312 

.0475793 

5 

-.7554044 

.0623145 

6 

-.6178762 

. 0747980 

7 

-.4580168 

. 0845783 

8 

-.2816036 

.0913017 

9 

-.0950125 

. 0947253 

10 

.0950125 

.0947253 

11 

.2816036 

.0913017 

12 

.4580168 

.0845783 

13 

.6178762 

. 0747980 

14 

. 7554044 

.0623145 

15 

.8656312 

.0475793 

16 

.9445750 

.0311268 

17 

. 9894009 

.0135762 1 


Pg order with an Sjg quadrature was adequate for obtaining the neutron dose rate in 
water. The results reported herein use an S^g Gauss- Legendre quadrature set 
(ref. 13) which has also been shown to be adequate for obtaining the fast neutron dose 
rate in a water medium. This quadrature set is listed in table III. 


ESTIMATION OF A VIRGIN SOURCE SPECTRUM 
Review of Small Source Spectra 

A review article by Hanson (ref. 14) covers much of the work on (a, n) sources up 
to about 1958. The dominant reaction in beryllium for alpha energies up to about 6 MeV 
is (refs. 15 to 17): 


9 Be + 4 He - 12 C + n + 5. 708 MeV 


1 9 

The C nucleus can be left in the ground state or in excited states at 4. 43 or 7. 65 MeV. 


10 


2 10 

For polonium 210 ( Po) alphas (E^ = 5. 30 MeV), the maximum neutron energy is about 
10. 9 MeV. 

Although measurement of the total 9 Be(cv, n) cross section (refs. 15 to 17) is reason- 
ably complete and the partial cross sections leading to the ground state and first excited 
12 

state of C have been measured (ref. 18), calculations of source spectra are not too 

satisfactory (ref. 19). Transitions to other than the ground and 4. 43-MeV states of C 

are accompanied by low energy neutrons so these partial cross sections are difficult to 

measure. St. Romain, Bonner, Bramblett, and Hanna (ref. 20) using Bonner spheres 

observed a group of neutrons at about 0. 3 MeV for alpha particles with about 5-MeV 

12 

energy. These neutrons were in addition to those from the 7. 65-MeV level of C which 
were not abundant (ref. 21). It was suggested (ref. 20) that low energy neutrons may be 
produced from endoergic breakup reactions in Be 

4 He + 9 Be - 4 He + [ 9 Be* - 8 Be + n - 1. 665 Mev] 


or 


[ 9 Be* - 2( 4 He) + n - 1. 564 Mev] 

9 

where the excited level in Be is at 2. 43 MeV or possibly at 1. 70 MeV. 

210 

In this work, reported spectra of several small Po-Be sources were studied to 
obtain what is thought to be a reasonable estimate of the virgin or uncollided neutron 
source spectrum. The virgin spectrum of neutrons from either plutonium 239-beryllium 
(23 9 p u _Be) (E^ = 5. 15 MeV) or ^^Am-Be (E^ = 5. 48 MeV) should be very similar to 
^* 9 Po-Be (E^ = 5. 30 MeV). However, the number of low energy neutrons from the en- 
doergic fraction increases with alpha energy. So the fraction of low energy neutrons in 

241 

a virgin source may be expected to be greatest in Am-Be sources. 

2 10 

The spectra of neutrons from Po-Be sources have been measured many times 
by a variety of methods. An early measurement is that of Whitmore and Baker (ref. 22) 
who used photographic emulsions. Notarrigo, Parisi, Ricamo, and Rubbino (ref. 23) 
used emulsions and reviewed the measurements up to 1961. They also made a calcula- 
tion of the source spectrum on an isotropic (a, n) angular cross section hypothesis. 

Haag and Fuchs (ref. 24) in 1962 used a stilbene crystal spectrometer following the me- 
thod of Broek and Anderson (ref. 25). One of the latest measurements, by Freestone, 
Muckenthaler, Straker, and Henry (ref. 26) at Oak Ridge (1967), used a similar tech- 
nique but with the NE-213 liquid scintillator (ref. 27). They show essentially identical 

23Q 210 241 

results for spectra of small Pu-Be, Po-Be, and Am-Be sources. 

Each method of measurement has sources of error. In the emulsion method, emul- 
sion shrinkage upon development, loss of track length through the surface, and difficulty 
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in identifying and measuring extremely short low energy track lengths are typical prob- 
lems. In a stilbene scintillator, anisotropy of light output with track orientation occurs. 
In an NE-213 scintillator, self-absorption of light and inability to distinguish between 
neutron and gamma induced events at small light output are problems. Both stilbene and 
NE-213 exhibit a strong nonlinearity of proton light response with energy. The shape of 
Haag and Fuchs’ spectrum at low energy may reflect an over compensation for this ef- 
fect. None of the methods are satisfactory much below 1 MeV. 

Figure 4 shows several measured Po-Be neutron spectra normalized to a yield of 

0. 88 neutron over the energy range from 1 to 11 MeV; a yield of 0. 12 below 1 MeV will 

be shown to be necessary for the interpretation of the present source experiments. All 

show broad maxima from 3 to 5 MeV. Most of these neutrons come from reactions 
12 

leaving C in the 4. 43 -MeV excited state. Smaller maxima near 7 and 9. 5 MeV arise 
out of anisotropy of the ground state branch (ref. 17). Freestone, Muckenthaler, 
Straker, and Henry (ref. 26) and Whitmore and Baker (ref. 22) are in best agreement, 
showing small differences. The spectrum of Haag and Fuchs is somewhat harder than 
the others. 

Various authors have suggested the presence of a large contribution to the total 
source spectrum below 1 MeV. St. Romain, Bonner, Bramblett, and Hanna studied a 
Pu-Be source that contained 13 grams Pu and probably 7 grams Be, By use of the 
Bonner sphere response functions, they assigned about 16 percent of the neutrons below 
1 MeV in this source, with a peak at 0. 3 MeV. Contributions to this portion of the spec- 
trum from scattering and (n, 2n) reactions in beryllium are estimated to reduce the frac- 
tion of low energy virgin neutrons in this source by about 3 percent of the total. Geiger 
and Hargrove (ref. 28) measured the spectrum of a small Am -Be source that contained 
3 grams of Be alloyed with 0. 23 grams of Am, enclosed in stainless steel 0. 13 centime- 
ter thick. They calculated 17 percent of the neutrons in their spectrum below 1. 5 MeV. 
These data indicate that a virgin Po-Be spectrum contains a significant fraction below 
1 MeV. 

A preliminary virgin spectrum was obtained from the data of figure 4. The shape 
from 1 to 11 MeV was constructed from the spectra shown by weighting the Whitmore 
and Baker and Freestone spectra twice as much as the other two measured spectra. The 
fraction from 0 to 1 MeV followed the shape given by St. Romain from 0 to 0. 6 MeV. 

This is shown in figure 4 as a dashed segment from 0 to 1 MeV. These yield shapes 
above and below 1 MeV were used to obtain a virgin spectrum that was compatible with 
the present experiments for the Am-Be source in air and in water. 


Derivation of a Consistent Virgin Source Spectrum 

An uncollided or virgin spectrum was constructed for the Am-Be source by means 
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of S n transport calculations for two experiments using the preliminary virgin spectrum 
of the previous section as a starting point. In the first experiment a leakage or collided 
spectrum for the source in air was measured using an NE-213 proton recoil scintillator. 
In the second experiment the indium resonance energy flux distribution was measured for 
the source in an effectively infinite water medium. By an iterative procedure using the 
S n method, a virgin spectrum was obtained which satisfied both experiments. This vir- 
gin spectrum was constrained not to exceed the estimated limits of error of published 
small source spectra measured by several independent techniques. 

Source in air. - The leakage spectrum for the Am-Be source in air was measured 
using an NE-213 proton recoil scintillation spectrometer. In the experiment, the Am-Be 
source was mounted 370 centimeters above an earthen mound to minimize wall and 
ground scattering. The top of the mound is about 9 meters long by 6 meters across and 
about 6 meters above ground. The source was supported by a 0. 64-centimeter-diameter 
brass rod positioned in an aluminum tube. 

The measured Am-Be source leakage spectrum down to 0. 5 MeV, obtained with an 
NE-213 proton recoil scintillation spectrometer and analyzed by the FERDOR computer 
program (ref. 29), is shown in figure 5. The solid histogram in figure 5 is the calcu- 
lated leakage spectrum using the yield shapes above and below 1 MeV as shown in fig- 
ure 4 and Perkins' evaluation for the (n, 2n) cross sections for the source beryllium. 

The required virgin spectrum is shown as the dashed line. 

The agreement of the calculation with the measurement of the flux per unit energy 
and per source neutron at 2 meters from the Am-Be source in air is very good over the 
energy range from 2 to 8 MeV. For the small neutron yields from 8 to 11 MeV, the 
calculation is lower than the measured data. Below 2 MeV the calculation is in good 
agreement with the measured data. Thus, the virgin source spectrum listed in table IV 
along with Perkins' Be(n, 2n) data gives a calculated leakage spectrum using the S n me- 
thod which fits the measured data quite well. 

Source in water . - Slowing down of neutrons in water is well understood. The spa- 
tial distribution of the indium resonance energy flux can provide additional evidence of 
the validity of the assigned source energy spectrum and the calculational method. 

Source neutrons of a given energy give a characteristic spatial distribution of the indium 
resonance energy flux. This can be seen from the distributions plotted in figure 6 for a 
number of point monoenergetic sources in water. The energy of the various monoener- 
getic sources corresponds to the average energy of a number of the energy groups listed 
in table II. The flux age for the various distributions are also listed in figure 6. The 
calculations were performed using PgS^g multigroup transport theory and correspond 
to a source strength of 1 neutron per second. 

Indium resonance energy flux distributions calculated for point monoenergetic 
sources having the average energy of each group listed in table II, result in an indium 
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TABLE IV. - FRACTIONAL ESTIMATED VIRGIN 


SOURCE AND GROUP AVERAGED Be(n, 2n) 
CROSS SECTIONS FOR 30-GROUP SPLIT 


OF TABLE II 


Total Be(n, 2n) cross sec- 
tion, b 


GAM-II data 


Perkins 1 data 


Group 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 


0. 559 
. 569 
. 579 
. 592 
. 600 
.605 
.606 
.603 
.600 
. 590 
. 565 
.376 
.026 
.00001 

0 


0. 528 
. 583 
. 594 
. 596 
. 596 
. 600 
. 597 
. 597 
. 593 
. 582 
. 533 
.296 
.039 
.0007 

0 


Virgin source 
fraction 


0.00862 
.02222 
.03534 
. 05469 
. 05650 
.06156 
.07081 
.07632 
.07572 
.07390 
. 12078 
. 09703 
.07858 
. 02786 
. 04087 
.0469 
.0433 
. 006 
. 003 

0 


resonance energy flux distribution in water <p( r) for a distributed source spectrum 
given by 


V(r) V> g (r) 

g 


( 3 ) 
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In equation (3) the quantity f is the source fraction in group g and (r) is the indium 

§ o 

resonance energy flux at distance r from a point source having the average energy of 
group g. 

Although the spatial distributions shown in figure 6 can be used to obtain indium 
resonance energy flux distributions in water for a given source spectrum, direct multi- 
group S n calculations are made that account for the source size and composition as 
well as the stainless steel container. An indium resonance energy flux distribution from 
a 30-group PgS^ transport calculation using the virgin source spectrum listed in ta- 
ble IV is shown in figure 7. Figure 7 also shows the normalized count rates multiplied 

o 

by the square of the radius C 0 r for the Am-Be source in water. 

Both the measured and calculated indium resonance energy flux distributions shown 
in figure 7 are normalized such that the volume integral of the slowing down density is 
unity. The total number of neutrons Q slowing down past indium resonance energy is 
given by 

/ oo 

q(r)dV (4a) 

where the slowing down density q(r) equals the product of the slowing down power of the 
medium times the flux. The quantity dV is the element of volume. For the experi- 
ments the measured count rate C is taken as proportional to the flux so that equation 
(4a) may be rewritten as 



where k(r) is the ratio of the slowing down power of the medium at r to that of water. 

In table I are given the normalized experimental saturated front and back averaged 
count rates of the indium foils C_ measured for the Am-Be source in water. These 

S 

data have also been corrected for foil size effect according to equation (2). The nor- 
malization factor n = 1. 837x10” is the factor required to make Q' = 1 in the experi- 
ment with the source in water. This same normalization factor is used in the treatment 
of the other experiments in water with beryllium shells around the source. As indicated, 
the calculations are normalized in a similar manner. 

2 

The logarithm of the normalized count rates Cr data shown in figure 7 can be 

S 

fitted by a straight line within the experimental error at distances greater than 25 centi- 
meters. If thermal neutron activation data for indium, not listed in table I, are normal- 
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ized between 45 and 65 centimeters to the indium resonance energy activity, the straight 
line fit extends out to 100 centimeters from the source. The measured data are thus fitted 
by the approximation 

(5) 

By using a least squares technique, the characteristic relaxation length X obtained is 
9. 54±0. 10 centimeters. This least squares fit to the measured data beyond 25 centime- 
ters obtained from equation (5) is shown as the dashed line in figure 7. 

In figure 7 the solid line shows the flux distribution from the 30-group PgS.^. trans- 
port calculation. The virgin source spectrum used was obtained by first using a virgin 
spectrum consistent with the leakage spectrum for the source in air and then modifying 
it to obtain agreement with the experiment by means of the PgSjg calculation of the 
flux distribution in water. As mentioned before, the source fraction in each group was 
subject to the constraint of being within the estimated bounds of the published small- 
source measurements. The virgin source fraction required below 1 MeV for best fit 
was almost exactly 12 percent, similar to the original estimate for the source in air. 

Calculated and measured parameters such as the flux age and higher spatial mo- 
ments can be compared for a finite sized source. The m th even spatial moment r m 
for the indium resonance energy flux distribution cp( r) in an infinite spherical medium 
is defined as 


C r = A exp 
s 


(i 


— ^ r m <p{r)dV 

r = — 

/ CO 

(p(r)dV 


(m = 0, 2, 4, 6, . . .) 


( 6 ) 


2 

The flux age r is one- sixth of the mean square of the distance r from the source at 
which neutrons possess indium resonance energy. In table V are some of the calculated 
and measured parameters for the source in water. These parameters also indicate the 
closeness of the fit established for the source in water with calculations based on the 
virgin source spectra. For example, the calculated flux age j is 55. 5 square centi- 
meters as compared to the measured 56. 0±1. 2 square centimeters for the finite sized 
Am-Be source in water. 

Calculated neutron flux spectra at a number of radii of figure 7 are shown in figure 8 
for the source in water. The mean energy of the flux spectrum at each radius is shown 
by the symbol on the curve. The mean energy and flux spectrum change only slowly with 
radius past about 10 centimeters. 


16 



TABLE V. - PARAMETERS OF EXPERIMENTAL 


AND CALCULATED INDIUM RESONANCE 
ENERGY FLUX DISTRIBUTIONS FOR 
Am-Be SOURCE IN WATER 


1 

i 

Experimental 

Calculated 

Second moment, 
r 2 xl0" 2 , cm 2 

3. 360 ( a l. 1) 

3. 331 ( b 3. 363) 

Fourth moment, 
r 4 xl0' 5 , cm 4 

3. 92 ( a 7. 6) 

3. 93 ( b 4. 02) 

Sixth moment, 
r 6 xl(T 8 , cm 6 

10. 81 ( b 25. 9) 

11. 23 ( b ll. 63) 

Relaxation length, 
A, cm 

9. 54t0. 10 


Flux age, r, cm^ 

56. 0±1. 2 

55. 5 ( b 56. 1) 


a Percent extrapolation for numerator (percent ex- 


trapolation for denominator ~0. 1). 
kpgSgg calculation. 

Accuracy of Calculation! Procedure 

A few remarks are in order concerning the validity of the comparison of the exper- 
imental leakage spectrum in air with the measured flux distribution in water. This com- 
parison rests upon the assumption that the S n calculation represents the spatial flux 
distribution correctly. While the S n calculation can do this, the fineness and type of 
the angular quadrature set, the scattering order of the cross sections, the number of 
energy groups, and the method of averaging the cross section data are important when 
detailed agreement is sought. The 30-group P 3 S 16 calculation used in this report is 
based on the limited storage capability of the IBM-7094 computer. It was possible, for 
the bare source in water, to perform a P 3 S 32 calculation. This has the effect of rais- 
ing calculated water fluxes at large distances from the source. Increasing the spatial 
mesh interval, beyond 10 centimeters, from 1 to 1. 5 centimeters caused a lowering of 
the calculated fluxes at distances past 20 centimeters of more than 2 percent. Finer 
spatial mesh interval and greater S n angular quadrature scattering order thus tend to 
increase the calculated flux at large distances for a given assumed source spectrum. 

For comparison, removal of 1 percent of the total source neutrons from the group 
above 10 MeV will lower the calculated resonance flux beyond 70 centimeters about as 
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much as the difference between the solid line in figure 7 and the dashed line. It is be- 
lieved that there is not a significantly larger fraction of neutrons in the first three 
groups than estimated in the solid line of figure 5. 

The measured water distribution is not very sensitive, however, to changes in the 
source that merely shift small fractions to neighboring groups. The sensitivity of a 
change in the fraction below 1 MeV is of some interest, since this fraction was not in- 
dependently determined. Removal of 2 percent of the total source from group 16, cen- 
tered at 0. 6 MeV, and increase of 2 percent in group 8 at about 5. 2 MeV will adequately 
represent such a shift. The calculated distribution in figure 7 decreases by a maximum 
of 2. 9 percent near a radius of 8 centimeters, is unchanged near 14 centimeters, and 
increases by as much as 3. 6 percent from 30 to 40 centimeters. The uncertainty in the 
measured points is somewhat less in the respective ranges. The estimated virgin 
source fraction below 1 MeV, using the Perkins' evaluated (n, 2n) data is then 12. 0±2. 0 
percent. 

COMPARISON OF CALCULATED RESULTS USING TWO EVALUATED Be(n,2n) 
DATA SETS WITH MEASUREMENTS IN WATER FOR SOURCE 
ENCLOSED WITHIN BERYLLIUM SPHERES 

Measurement of the flux distribution in water at indium resonance energy for the 
Am-Be neutron source enclosed by a beryllium shell yields several quantities which can 
be calculated. The spatial neutron flux distribution and the number of neutrons slowing 
down are the most important. Spatial moments of the flux distribution can also be cal- 
culated. Two shells of beryllium were used in the experiments. The smaller shell was 
1. 91 centimeters thick and had an outside diameter of 10. 16 centimeters. The larger 
shell was 7. 0 centimeters thick and had an outside diameter of 20. 32 centimeters. 


Detailed Shapes of the Indium Resonance Energy Flux Distributions 

Figures 9(a) and (b) show two experimental flux distributions in water along with 
their corresponding multigroup transport calculations. The S n calculations in each 
case use an S^g angular quadrature order, Pg elastic scattering order, Pq (n, 2n) scat- 
tering order for Be, and 30 energy groups. Two S n calculations were performed for 
each experiment. One of the calculations used Perkins’ Be(n, 2n) evaluation to compute 
Pq (n,2n) down scattering transfer cross sections with the AGN-SIGMA computer pro- 
gram. The other calculation used Be(n, 2n) Pq down scattering transfer cross sections 
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from the GAM-II computer program data files. Both calculations used the virgin source 
spectrum listed in table IV. 

Differences between the calculations with Perkins’ (n, 2n) data for Be and GAM-II 
Be(n, 2n) data are significantly large for the 7. O-centimeter-thick Be shell flux measure- 
ment shown in figure 9(a). The flux measurement for the 1. 91-centimeter-thick Be 
shell shown in figure 9(b) indicates similar differences for the two calculations. In the 
7. O-centimeter thick Be shell the calculation which used GAM-II Be(n, 2n) data does not 
fit the measurements. Only from about 10 to 15 centimeters is there agreement. From 
15 centimeters on this calculated distribution is much larger in magnitude than the mea- 
sured distribution. For this 7. O-centimeter-thick Be shell the calculation which used 
Perkins' Be(n, 2n) data does fit the measurements. From 15 to 30 centimeters this cal- 
culated distribution varies by only a few percent from the measured distribution. From 
30 centimeters to the end of the measured data the calculated distribution is slightly 
lower than the measured distribution. The calculation is higher than the measurement 
from about 10 to 15 centimeters. Perkins' Be(n, 2n) evaluation includes alternative 
models for the multibody component. The calculations shown here as using Perkins' 
Be(n, 2n) data use the three-body phase space model for the multibody component. A 
calculation using Perkins' Be(n, 2n) data which uses the cluster model for the multibody 
component gives a generally poorer fit to the measured data for the 7. O-centimeter- 
thick Be shell experiment. The calculation with the cluster model component is lower 
than the calculation with the three-body component between 15 and 35 centimeters, the 
maximum difference being 6 percent at about 20 centimeters. 

Remarks that were made previously concerning the detailed agreement of the trans- 
port calculation with the measurement for the bare source in water hold also for the 
calculations with the beryllium shell enclosed source in water. For example, to the 
extent that the calculated distribution is higher relative to the measurement in water at 
large distances from the source, the same should be true for the beryllium shell en- 
closed source in water results. In addition, the effect of neglecting the higher order 
Legendre components of the Be(n, 2n) down scattering transfer cross sections was es- 
timated by calculating the indium resonance energy flux distribution using the 2. 43- 
MeV excited state model for the total Be(n, 2n) cross section. Calculation of the flux 
distribution outside the 7. O-centimeter-thick Be shell, but near the shell, was about 
2 percent lower when the P^ down scattering transfer cross sections for the Be(n, 2n) 
reaction were included, which accounts for about 20 percent of the difference between 
the calculated and measured values in this region. 

Although neutrons of a given energy slow down throughout the entire volume, the 
higher energy neutrons are more important at larger radii from the source. The cal- 
culation using the GAM-II Be(n, 2n) data shows more neutrons slowing down at radii 
greater than 15 centimeters than either the experiment or the calculation using the 
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TABLE VI. - PARAMETERS OF EXPERIMENTAL AND CALCULATED INDIUM RESONANCE ENERGY 


FLUX DISTRIBUTIONS FOR Am-Be SOURCE ENCLOSED BY Be SHELL FOLLOWED BY WATER 



1.9- 

cm-thick Be shell 

7.0- 

cm-thick Be shell 


Be(n, 2n) data 

Experimental 

Calculated 

Be(n, 2n) data 

Experimental 

Calculated 

Second moment, 

Perkins’ 

3. 069 ( b l. 4) 

3. 074 

Perkins’ 

2.419 ( b l. 4) 

2. 397 

r 2 xl0' 2 , cm 2 

GAM-n 

3. 092 ( b l. 4) 

3.211 

GAM-II 

2. 550 ( b l. 4) 

2.696 

Fourth moment, 

Perkins' 

3. 232 ( b 9. 5) 

3. 391 

Perkins’ 

2.003 ( b 10. 6) 

1.996 

r"S<10~^, cm^ 

GAM-n 

3. 257 ( b 9. 5) 

3. 537 

GAM-n 

2. 140 ( b 10. 6) 

2. 322 

Sixth moment, 

Perkins' 

8. 058 ( b 29. 6) 

9.273 

Perkins’ 

4. 625 ( b 33. 4) 

4. 917 

r^xlCT^, cm^ 

GAM-II 

8. 121 ( b 29. 6) 

9. 577 

GAM-II 

4. 947 ( b 33. 4) 

5. 582 

Relaxation length, 


8. 96 ±0. 10 



9. OOtO. 10 


X, cm 








a Two experimental values for each moment correspond to use of calculated flux distribution in Be shell 


based on indicated Be(n, 2n) data set. 

^Percent extrapolation for numerator (percent extrapolation for denominator ~0. 1). 


Perkins’ Be(n, 2n) data. This emphasizes the fact that the Be(n, 2n) Pq down scatter- 
ing transfer matrix used by the GAM-II evaluation retains more high energy neutrons 
than does Perkins' evaluation. 

The experimental moments for these cases are given in table VI. Moments obtained 
from the calculated distributions for the two Be(n,2n) evaluations are also shown in 
table VI. As expected, agreement with experimental moments is better using the calcu- 
lation with the Perkins' Be(n, 2n) data than with the GAM-II Be(n, 2n) data. 

Calculated neutron flux energy spectra for the 7. O-centimeter-thick Be shell exper- 
iment are shown in figure 10. The calculation based on Perkins' Be(n, 2n) data is shown 
by the solid lines, while the calculation with GAM-II Be(n, 2n) data is shown by the dashed 
lines. The average energy is also given for each spectrum shown. 


Enhancement Due to Beryllium 

The volume integral of the indium resonance energy slowing down density is a meas- 
ure of the number of neutrons slowing down past 1. 45 eV per unit time per virgin source 
neutron. In addition to the production of neutrons by the (n, 2n) reaction in Be, there are 
high energy neutron absorptions in every isotope present except hydrogen. The net num- 
ber of neutrons at 1. 45 eV for the two Be shell experiments are related to a quantity 
called enhancement. The enhancement is defined here as the fractional increase in the 
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TABLE VII. - NEUTRON ENHANCEMENT FOR Am-Be SOURCE ENCLOSED BY Be SHELL FOLLOWED BY WATER 


Experiment Be data set Source Number of neutrons absorbed above Net number of neutrons Relative enhancement, percent 

description augmenta- 1. 45 eV slowing down past T” 

a I ,, i e y Calculated Experimental 0 

Be in Source con- Be shell Water 
source tainer and 


spacer if 
present 


Source in water 

Perkins’ 

1.0486 

0.0064 

0.0007 


0.0266 

1.0149 

_ 



GAM-II 

1. 0472 

.0065 

.0007 


.0266 

1.0134 

— 


Source with 1. 9- 

Perkins' 

1. 1572 

0. 0067 

0.0011 

0.0166 

0.0220 

1.1108 ( c 0. 0089) 

9.4 

7. 3±2. 0 

cm-thick Be shell 

GAM-II 

1. 1587 

.0069 

.0010 

. 0183 

.0225 

1. 1100 ( c 0. 0080) 

9. 5 

7. 2-k2. 0 

Source with 7.0- 

Perkins’ 

1. 3578 

0. 0068 

0.0014 

0. 0540 

0.0143 

1.2813 ( c 0. 0816) 

26. 3 

23. 5t2. 5 

cm-thick Be shell 

GAM-II 

1.3817 

.0071 

.0012 

.0651 

.0150 

1.2933 ( c 0. 0675) 

27.6 

22. U2. 5 


a Source augmentation is total virgin source neutron (1. 0) plus number of neutrons produced by (n, 2n) reaction in Be. 

^For 1. 9-cm-thick Be shell experiment, estimated measured enhancement is taken to be 7. 3±2. 0 percent; for 7. 0-cm-thick Be shell experi- 


ment, estimated measured enhancement is taken to be 22. 8±2. 5 percent. 
c Net number of neutrons slowing down past 1. 45 eV in Be shell. 


number of neutrons slowing down past indium resonance energy for the source enclosed 

in a Be shell in water as compared to the source alone in water. Table VII shows the 

calculated neutron productions, absorptions, the net number of neutrons slowing down 

past 1. 45 eV, and the enhancements. The fission and absorption cross sections of 
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Am are ignored in these calculations. The energy of 1. 45 eV corresponds to the high 
energy boundary of group 28 listed in table II. 

It should be noted that the indium resonance energy flux distribution was not meas- 
ured in the Be shells. Therefore, to obtain the experimental enhancement, the calcu- 
lated flux distributions were used for the Be shells. Although the contribution to the flux 
volume integral from the Be shells and source are large fractions of the total flux volume 
integral at indium resonance energy, their contribution to the slowing down integral is 
much less amounting to about 6 percent in the 7. O-centimeter Be shell experiment. The 
reason for this is that the ratio of the slowing down power of Be to that of water is 0. 106. 
This value was based on the elastic scattering cross sections for Be and water. Equa- 
tion (4b) is then used to calculate the total slowing down density at indium resonance 
energy. 

If the enhancements in table VII are considered, not much significance can be at- 
tached to the fact that there is a slight difference between the two calculations of the en- 
hancement in the experiment with the 1. 9-centimeter-thick Be shell. In the experiment 
with the 7. O-centimeter-thick Be shell, however, the calculated enhancements differ: 

26. 3 percent using Perkins' Be(n, 2n) data and 27. 6 percent using GAM-II Be data. When 
the experimental enhancements are considered, two values are shown because of the 
need to use calculated flux distributions in the Be shell. A reversal can be seen between 
the experimental enhancements and the calculated enhancements (i. e. , 23. 5±2. 5 and 
22. 1±2. 5 percent) for the experimental values using calculated Be shell flux distributions 
based on Perkins’ and GAM-II Be(n, 2n) data, respectively. This results because there 
was, of course, only one experiment and the neutron flux distribution as calculated using 
Perkins’ Be(n, 2n) data is higher than that using the corresponding GAM-II data. 

The errors quoted on the experimental enhancements are the estimated standard 
error of the measurements and do not reflect the uncertainty in the flux in the Be shell 
and the source. Thus, in the 1. 9-centimeter Be shell experiment the calculated value 
is just outside the 7. 3±2. 0 percent estimated for the experimental value. In the 7. 0- 
centimeter-thick Be shell experiment the enhancement calculated using Perkins’ Be(n, 2n) 
data is just outside the error limit for the estimated experimental value. However, the 
calculated enhancement using GAM-II Be(n, 2n) data considerably exceeds the error limit 
for the experimental value. Thus, calculations using Perkins' Be(n, 2n) data are in bet- 
ter agreement with the experimental measurements than calculations using GAM-II 
Be(n, 2n) data. However, the results for these transmission experiments with Be shells 
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suggests the need for some revision of Perkins’ evaluation for the (n,2n) reaction in Be. 


CONCLUSIONS 

A number of conclusions concerning two Be(n, 2n) cross section evaluations result 
from this study of the measured indium resonance energy flux distributions in water due 
to neutrons transmitted through thick Be shells enclosing an americium-beryllium 
(Am -Be) source: 

1. Multigroup P 3 S 16 transport calculations with a virgin source spectrum and 
Be(n, 2n) cross section data from Perkins’ evaluation gave good agreement with indium 
resonance energy flux distribution measurements in water for the Am-Be source en- 
closed by Be shells. 

2. Calculations with Be(n, 2n) cross section data from the GAM-II evaluation gave 
poor agreement with indium resonance energy flux distribution measurements in water 
for the Am-Be source within the Be shells. 

3. For the Am-Be source in the 1. 91-centimeter-thick Be shell, the computed en- 
hancements are 9. 4 and 9. 5 percent for calculations using the evaluation due to Perkins 
for the Be(n, 2n) cross sections and the evaluation used by the GAM-II program, respec- 
tively. The experimental enhancement is 7. 3±2. 0 percent. 

4. For the Am-Be source in the 7. 0- centimeter -thick Be shell, the computed en- 
hancements are 26. 3 and 27. 6 percent for calculations using the evaluation due to Per- 
kins for the Be(n, 2n) cross sections and the evaluation used by the GAM-II program, 
respectively. The experimental enhancement is estimated to be 22. 8±2. 5 percent. The 
calculation using Perkins* evaluation for the Be(n, 2n) cross sections is better than that 
using the GAM-II evaluation. 

5. These results show that the calculation of the spatial distribution of inditun re- 
sonance energy neutrons and of the enhancement, as compared to the measurements, 

is significantly better using Perkins’ evaluation for the (n, 2n) data for Be than the eval- 
uation used by the GAM-II program. Since current ENDF/B cross sections are based on 
the same evaluation as used to formulate the GAM-II cross sections, it is suggested that 
the Perkins’ evaluation be used to update cross sections for beryllium. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 26, 1970, 

129-02. 
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Figure 8. - Neutron spectra for Am “Be source in water medium. 
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(a) Am-Be neutron source enclosed by 7.0-centimeter-thick Be shell followed by water. 
Figure 9. - Flux distribution at indium resonance energy times radius squared. 
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Figure 10. - Neutron spectra for Am-Be source enclosed in 7.0-centimeter-thick Be shell and followed by water. 
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